Role of spin-orbit coupling on the electronic structure and properties of SrPtAs 
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The effect of spin-orbit coupling on the electronic structure of the layered iron-free pnictide su- 
perconductor, SrPtAs, has been studied using the full potential linearized augmented plane wave 
method. The anisotropy in Fermi velocity, conductivity and plasma frequency stemming from the 
layered structure are found to be enhanced by spin-orbit coupling. The relationship between spin- 
orbit interaction and the lack of two-dimensional inversion in the PtAs layers is analyzed within 
a tight-binding Hamiltonian based on the first-principles calculations. Finally, the band structure 
suggests that electron doping could increase T c . 
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I. INTRODUCTION 

The recent discovery of superconductivity in pnictides 
has attracted extensive attention owing to their surpris- 
ingly high T c }- While the highest T c so far is 56 K for 
GdFeAsO^ a consensus on the pairing mechanism has 
not yet been reached £ This class of materials share a 
common crystal structure, that is, the Fe square lat- 
tice. While most of the superconducting pnictides are 
Fe-based, pnictides without iron also exhibit supercon- 
ductivity, although T c is drastically lower than with iron. 
Recently, another superconducting pnictide, SrPtAs, has 
been discovered, which is the first non-Fe based super- 
conductor with a hexagonal lattice rather than square 
latticed Although T C —2AK is lower than those of Fe- 
based pnictides, it possesses interesting physics associ- 
ated with its hexagonal crystal structure. 

SrPtAs crystallizes in a hexagonal lattice of ZrBeSi 
type with space group PG^/mmc (No. 194, — the 
same (non-symmorphic) space group as the hep structure 
— with two formula units per primitive cell. As depicted 
in Fig. [IJa) , its structure resembles^ that of MgB2 (with 
the symmorphic space group P6/mmm, No. 191, Dqu), 
with a double unit cell along the c axis: the boron layers 
of MgB 2 are replaced by PtAs layers, rotated by 60° in 
successive layers (responsible for the non-symmorphic na- 
ture) and Mg is replaced by Sr. Although the crystal as a 
whole has inversion symmetry (as do the Sr atoms), the 
individual PtAs layers lack two-dimensional inversion. 

Thus, SrPtAs differs from MgB2 in two significant 
ways: (i) it exhibits strong spin-orbit coupling (SOC) 
at the Pt ions and (it) the PtAs layers individually 
break inversion symmetry, exhibiting only C 3 „ (or D^h 
if ^-reflection is included) symmetry. These two prop- 
erties play an important role in determining the band 
structure and also affect the superconducting state. As- 
suming that the superconductivity is largely determined 
by the two-dimensional PtAs layers, the lack of in- 
version symmetry in the individual layers (which we 
call "broken local inversion symmetry") opens up the 



possibility to see the unusual physics associated with 
non-centrosymmetric superconductors.^ With large spin- 
orbit coupling, nominally s-wave non-centrosymmetric 
superconductors exhibit spin-singlet and spin-triplet 
mixing)^ enhanced spin susceptibilities^^ enhanced 
Pauli limiting fields^ non-trivial magnetoelectric effects 
and Fulde-Ferrell-Larkin-Ovchinnikov(FFLO)-like states 
in magnetic fields^— and Majorana modest The local 
inversion symmetry breaking, together with a SOC that 
is comparable to the c-axis coupling, suggests that SrP- 
tAs will provide an ideal model system to explore related 
effects in centrosymmetric superconductorsj 18 ' 19 

In this paper we discuss the electronic structure of SrP- 
tAs, including spin-orbit coupling, which was neglected 
in a previous theoretical study^ In Sec. UH we describe 
details of the calculations. The effects of SOC on the 
bands, the Fermi surface, density of states, and transport 
properties at the Fermi surface of SrPtAs are presented 
in Sec. IIII1 along with a tight-binding analysis. Finally, 
we suggest an enhanced T c might be possible via doping. 
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FIG. 1. (a) (Color online) Crystal structure of SrPtAs, where 
red, blue, and grey spheres denote Pt, As, and Sr atoms, 
respectively, (b) Brillouin zone of SrPtAs and high symmetry 
k points. 
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II. METHOD 

First-principles calculations are performed using 
the full-potential linearized augmented plane wave 
(FLAPW) method 21 ' 2 - 2 and the local density approxi- 
mation (LDA) for the exchange-correlation functional of 
Hedin and Lundqvisti^ Then, SOC is included by a 
second variational method^ Experimental lattice con- 
stants, a= 4.244A and c= 8.989A, are used^ Cutoffs 
used for wave function and potential representations are 
196 eV and 1360 eV, respectively. Muffin-tin radii are 
2.6, 2.4, and 2.1 a.u for Sr, Pt, and As, respectively. 
Semicore electrons such Sr Ap and As 3d are treated as 
valence electrons, which are explicitly orthogonalized to 
the core states^ Brillouin zone summations were done 
with 90 k points in the Monkhorst-Pack scheme,— while 
the density of states are obtained by the tetrahedron 
method. 28 Although most of the calculations were per- 
formed using LDA, some results were also done with the 
GGA as well.— In order to calculate the Fermi velocity 
and plasma frequency for in-plane and out-of-plane con- 
tributions, eigenvalues from self-consistent calculations 
are fitted by a spline method over the whole Brillouin 
zone^Sr— 



III. RESULTS 

The band structure of SrPtAs is presented in Fig. m 
along the symmetry lines shown in Fig. [ljb). Plots in 
the left (right) column are without (with) spin-orbit cou- 
pling, and plots in the upper and lower rows are the same 
but highlighted for As p and Pt d orbitals, respectively. 
Our energy bands and Fermi surfaces without SOC agree 
with those of Ref. [2(| ■ The main band of Sr 5s origin is 
located far above the Fermi level (Ep), consistent with 
Zintl's scheme^ that Sr donates electrons to the PtAs 
layer and behaves almost like an inert Sr 2+ ion. Without 
SOC, bands on the zone boundary face, the k z = tt/c 
plane(i/ — A — L — H), exhibit four fold degeneracy - 
two from spin and the other two from two different PtAs 
layers — as a consequence of the non-symmorphic trans- 
lations along the c-axis, just as for the hep structure. 
For the k z = plane (K — T — M — K), there is no such 
symmetry-dictated degeneracy due to the two equivalent 
layers, but instead the magnitude of splitting is propor- 
tional to the inter-layer coupling. With SOC, the bands 
change markedly: The four-fold degeneracy on the zone 
boundary face is reduced to a two fold pseudospin degen- 
eracy due to inversion symmetry, whereas bands along 
the A — L line keeps the fourfold degeneracy as a conse- 
quence of time-reversal symmetry^ 

In a simple atomic picture, the SOC Hamiltonian is 
H SO c — SL-cr, where S represents the strength of the SOC. 
Values of the SOC strength derived from the calculations 
are 6pt=0.32 eV for the Pt d orbitals and <5as=0.23 eV 
for As p orbitals, which are used in later discussions. The 
SOC splitting at the A point for Pt d XZl d yz (d xy , d x 2_ y 2) 



orbitals is 0.59 (0.54) eV, whereas for Asp XtV the splitting 
0.28 eV. 

To gain insight into the effect of spin-orbit coupling, we 
consider a simple tight-binding theory for a single PtAs 
layer. A single PtAs layer lacks a center of inversion 
symmetry and therefore allows an anti-symmetric spin- 
orbit coupling of the form 
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exists, where cj, ( c fc s ) creates (annihilates) an elec- 
tron with momentum k and pseudo-spin s, and <r de- 
note the Pauli matrices. Time-reversal symmetry im- 
poses gfc = — g fa- Here we find the form of g^ 
through a consideration of the coupling between the Pt 
d x 2_ y 2,d xy orbitals and the As p x ,p y orbitals. These 
orbitals give rise to the Fermi surfaces with cylindrical 
topology around the T-A line. On the Pt sites, we define 
states |d±, s >= \d x 2_ y 2, s > ±i\d xy , s > and on the As 
sites we define the states |p±, s >— \p x , s > ±i\p y , s > (s 
denotes spin). On a single Pt or As site, the spin-orbit 
coupling L ■ S has only L Z S Z with non-zero matrix ele- 
ments in these subspaces. This splits the local four-fold 
degeneracy into two pairs so that for Pt (As) \d+, f> and 
\d—,l> (|jH-,t> an d \p— ,4>) form one degenerate pair 
while |d+,|> and |d-,t> (\p—,t> an d \p+,l>) form 
the other degenerate pair. For both pairs, time-reversal 
symmetry is responsible for the degeneracy and we can 
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FIG. 2. (Color online) Band structure of SrPtAs (a),(c) with- 
out and (b),(d) with spin-orbit coupling. In (a) and (b), As 
p x ,y and p z orbitals are shown in red and green, respectively, 
and As s in blue. In (c) and (d), Pt {d xy ,d x 2_ y 2), (d xz ,d yz ), 
and d z 2 orbitals are presented in red, green, and blue, respec- 
tively. Contribution from Pt s is shown in purple. 
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label the two degenerate partners of each pair through a 
pseudo-spin index. We include a spin-independent near- 
est neighbor hopping between the As and Pt sites. This 



/ (d + S Pt 

e d - S Pt 

\ t(l + e lkT3 + e- !kT2 ) t(l + is*e lkT 3 + ue"^^) 



where v = e l2 "/ 3 , T ± = a(l,0), T 2 = a(-l/2, V3/2), 
T3 = a(— 1/2, — -\/3/2), i and i are the hopping param- 
eters between the As p and Pt <i orbitals, and Spt (6a s ) 
is the atomic SOC parameter for Pt d(As p) orbitals. To 
find an effective Hamiltonian for the Pt <i-orbitals, we as- 
sume that \ed — t P \ is the largest energy scale and treat 
all other parameters (t, i, 5pt, Sa s ) as perturbations. This 
yields the following effective Hamiltonian 

Up t = ^p t J k )Hp ttS (k)^p t ^(k), (4) 
k^s I 



Hpt,t — 



where h = 2[P + cos(27r/3)i 2 ]/(e rf - e p ), t 2 = 
-2sin(2^/3)t 2 /(e rf - e p ), and t 3 = -2tt/(e d - e p ), 
and the sum over i is over the three vectors Ti,T 2 , 
and T3. The Hamiltonian Hp t ^ yields the disper- 
sions e(k) = ii(k) ± y/(6 Pt + ^(k)) 2 + |t 3 (k)| 2 where 
ii(k) = h J2 t cos(k • T<), t a (k) = t 2 sin(k • T<), and 
i 3 (k) =t 3 [cos(k-Ti) + i/cos(k-T2) + i/*cos(k-T 3 )]. The 
SOC contribution that lifts the four-fold degeneracy into 
two bands i?k,± = ± |<?(k)| at k z = ir/c can be found 
to order 6 Pt /^S 2 pt + |t 2 (k)| 2 + |i 3 (k)| 2 and is given by 

/, n * t 2 5p t 

9[ >~ V^p t + l*2(k)| 2 + |t 3 (k)| 2 

x [sin(k ■ Ti) + sin(k • T 2 ) + sin(k ■ T 3 )] (6) 

where we have also included the contribution from 
Hpt,i(k). Note that this g leads to pseudo-spin in- 
teraction in Eq. [1] denoted by a z . We emphasize that 
this a z operates on pseudo-spin, not actual spin, the 
up and down pseudo-spin states are related by time- 
reversal symmetry (for example, the local states |d+,t> 



yields the following tight binding Hamiltonian in fc-space 
W = X>t(fe)tf.(k)*.(fe), (2) 

fe,s 

where * s (fe) = (c fe d+ s , c k ^ s , c fep+ s , c fc p _ s ) T , s is 
the spin label (s = F° r s P m u Pj Hf(k) is (to get 

Hi(k) change the sign of Sp t and 5a s ) 



t(l + ^*e- lkT3 + ^e lkT2 ) t(l + e- lkT3 + e* kT2 ) \ 
t(l + e- ikT i + e !k ' T2 ) t(l + iye- lkT3 + is*e lkT *) 
e p + 5 As 

e p - 5 As ) 

(3) 

I 

where \E , pt, s (fc) = (c^. d+ s ,Cf s , d _ s ) T , s is the spin label, 
and 



I 

and \d— ,4> form a pseudo-spin pair). The expres- 
sion for g clearly reveals how the interplay between the 
atomic SOC (<5pt) and the broken inversion symmetry 
(</(k) = — g(— k)) of a single PtAs layer leads to the rel- 
evant band SOC. Note that this single-layer band SOC 
will be of opposite signs for the two inequivalent PtAs 
layers in the unit cell. Further, the band spin-orbit split- 
ting will have additional contributions from terms of or- 
der 5pt/{€d — £p) and 5as/ {^d — £p) that were neglected in 
the above derivation. However, these additional contri- 
butions do not qualitatively change the results. Similar 
considerations apply for the other Pt <i-orbitals. 

The previous paragraph considered a single PtAs layer. 
To complete the description, the coupling between the 
two inequivalent layers must be included. For the Pt 
d x 2_ y 2,d xy orbitals considered above, the nearest neigh- 
bor inter-layer hopping matrix between \d+, s > states is 
(the same expression appears for \d— , s > states) 

e c (k) = t c cos(fc z c/2)(l + e~ lk - T3 + e ik - T2 ). (7) 

Including this inter-layer hopping leads to the following 



( e d + S Pt + h cos(k • + 1 2 £< sin(k • T<) i 3 [cos(k • Ti) + z/cos(k • T 2 ) + v* cos(k • T 3 )] \ 
^ < 3 [cos(k ■ Ti) + v* cos(k • T 2 ) + v cos(k • T 3 )] e d - S Pt + h J2 t cos(k • T<) - t 2 ^ sin(k • T,) J 
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Hamiltonian 

H±,s = ^±A k ){i e ±( k ) - A^oTu +fl(k) ■ (tt z 
k 

+Re[e c (k)}(ToT x + Im[e c (k)]aoT v 

where *±,,(fe) = (c ±k ^ s , c ±fe;i s , c ±k ^ 2 s , c ±ki2 s ) T , 
1,2 denote the two inequivalent PtAs layers, <ii (ji) are 
Pauli matrices that operate on the pseudo-spin (layer) 
space, e± = ti(k) ± ^/(i 2 (k)) 2 + |< 3 (k)| 2 , and t/(k) is 
given in Eq. ([6]) (the t z matrix describes the sign change 
of g on the two layers). This Hamiltonian can be di- 
agonalized with resulting dispersion relations e(k) = 
e± (k) ± ^/|e c (k)| 2 + g 2 (k) and each state is 2-fold degen- 
erate due to time-reversal symmetry (Kramers degener- 
acy) . Note that the tight binding theory described above 
suggests that eigenstates of S z are also eigenstates of the 
single electron Hamiltonian. However, inter-layer cou- 
pling terms can lead to additional terms that do not com- 
mute with S z . The band structure suggests that these 
terms are not large for the states near the Fermi surface. 

The SOC found in in Eq. ^ has opposite sign for the 
different layers as well as for the pseudo-spin direction. 
This is demonstrated in Fig. [3l where the band structure 
is resolved by layer and spin for H—L-H', where L is one 
of the time-reversal invariant momentum (TRIM) points, 
and the lines H-L and L-H' are related both by a mirror 
plane and, and more relevant to the discussion, by a com- 
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FIG. 3. (Color online) SOC splitting of bands, (a) and (b) 
[(c) and (d)] show components from Pt 5d orbitals for up- 
per [lower] Pt atoms. Left column[(a) and (c)] is for spin-up 
components and right[(b) and (d)] is for spin-down compo- 
nents. Radius of circles is proportional to the magnitude of 
the components. 



bination of inversion and reciprocal lattice vectors. In 
Fig. [3j spin up[(a)] and down[(b)] components of the Pt 
d orbital from the upper PtAs layer are marked in differ- 
ent colors, while components from the other PtAs layers 
are shown in (c) and (d). As expected from the tight- 
binding analysis, Figs. [Ha), (d) [and similarly for (b) and 
(c)] appear the same since they correspond to spatial in- 
version and time- reversal (opposite spin). Despite the 
presence of a global inversion center, the locally broken 
inversion symmetry in PtAs is evident in Fig. [3] where 
the spin degeneracy is broken in a single layer, for which 
the consequences are nontrivialJ^ (This result is not sur- 
prising or unexpected since in the limit that the coupling 
between PtAs layers vanishes, i.e., the separation goes to 
infinity, the single layer result must be recovered.) This 
spin separation, however, does not result in magnetism 
because of spin compensation in each layer. These anti- 
symmetric splittings (g(k) = — g(— k)) occur not only at 
the L point but also at the other TRIM points^ — T, A, 
M, and the A-L line — in the hep structure. 

As expected from the layer structure, anisotropy oc- 
curs in the Fermi surfaces, conductivity, and plasma fre- 
quency. Two-dimensional cross sections of the Fermi sur- 
faces are shown in Fig. 01 both without [(a)-(d)] and 
with [(e)-(h)] SOC. Contours in the k z = tt/c [(a),(e)] 
and k z = [(d), (h)] planes clearly exhibit the conse- 
quence of the crystal symmetry: hexagonal symmetry 
around T and A and trigonal symmetry around K and 
H. All sheets exhibit almost two-dimensional cylindri- 
cal features except for a small pocket around H. This 





FIG. 4. Cross section of the Fermi surface (a)-(d): without 
and (e)-(h): with spin-orbit coupling. Cross section along the 
zone boundary face, k z = n/c [(a),(e)] and at the zone center, 
k z = [(d), (h)]. (b),(c) are contours along the vertical plane 
shown with their corner points in the Brillouin zone, where 
(f),(g) are those with SOC included. Numbers next to Fermi 
surface sheets indicate band indices. 
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anisotropy is further exemplified by transport proper- 
ties which will be discussed later. Moreover, because 
of the symmetry-dictated degeneracy due to the non- 
symmorphic symmetry, there is no spin-orbit splitting 
along the time-reversal invariant direction, A-L, as seen 
in Fig. EJe). Comparing Figs.Htc) and (g), the SOC ap- 
pears to make the Fermi surfaces more cylindrical, conse- 
quently enhancing the two dimensional character of the 
Fermi surfaces. 

All the Fermi surfaces are hole-like after turning on 
SOC, in sharp contrast to other pnictide superconductors 
with two electron-like Fermi surfaces and two hole-like 
Fermi surfaces^ Instead of electron-like and hole-like 
Fermi surfaces, they are distinguished by orbital char- 
acter. Sheets around the T-A line consist of a orbitals 
of the PtAs layer, As p x>v and Pt d xy x 2- y z, while sheets 
around the K-H line are from ir orbitls, As p z and Pt 
d X z,yz- Two kinds of Fermi surfaces with different orbital 
character might give rise to a two energy gap supercon- 
ductor in SrPtAs. 

The anisotropy due to the layered structure is further 
manifested in the average Fermi velocities and plasma 
frequencies. Neglecting SOC, the in-plane and out-of- 
plane Fermi velocities are (v x ) 1 / 2 =3.72xl0 7 cm/s and 
(u2) 1 ^ 2 =l-02x 10 7 cm/s, respectively, and the plasma 
frequencies are D, XtV =5.70 eV and f2 z =1.57 eV. The 
anisotropy ratio, defined as a ratio of conductivities be- 
tween in-plane and out-of-plane components, is 13.3, as- 
suming an isotropic scattering rate. With SOC, the 
anisotropy is enhanced: (w 2 y ) 1 / 2 =3.76x 10 7 cm/s and 
(« 2 )V2 = 6.78xl0 6 cm/s; 0^=5.57 eV, n z =1.00 eV, and 
the anisotropy ratio increases to a much higher value 
of 30.8. The decrease of v z by 33% by SOC is consis- 
tent with the enhanced two dimensional character of the 
Fermi surfaces. Table U summarizes contributions from 
each Fermi surface. The largest contribution to the den- 
sity of states at the Fermi level, N(Q), comes from the 
34th band at around K, which is due to the low veloc- 
ity at the Fermi surface. The anisotropy ratio is usu- 
ally much larger than 1 except for the small hole pocket 
around H . 



TABLE I. Fermi surface properties with SOC included for 
each surface: Density of states, iV(0) (states/eV/spin); av- 
erage velocities, (v 2 ) 1 ^ 2 , (v 2 ) 1 ^ 2 (10 7 cm/s); anisotropy ratio, 
(v 2 )/{v 2 ); and plasma frequencies, Q x , Q, z (eV). The numbers 
31,32,33, and 34 represent band indices and F and K in the 
parenthesis represent the locations of the Fermi surface. 
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32 


33(r) 


34(r) 


33(K) 


34(K) 
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N(0) 


0.085 


0.107 


0.209 
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0.267 


0.943 
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(vl) 1/2 
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(v 2 ) 1/2 


1.03 


0.67 


0.47 


1.04 


1.37 


0.41 


0.678 


(vl)/(vl) 


46.9 
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FIG. 5. (Color online) Density of states (DOS) of SrPtAs: 
(a) total DOS, (b) s orbitals, (c) As p orbitals, and (d) Pt d 
orbitals. 



The total density of states (DOS) and orbital decom- 
posed partial DOS are presented in Fig. [S] The states 
around Ep arise mainly from As p and Pt d. The non- 
bonding Pt d z 2 bands are located at around -2.2 eV, and 
because they hybridize little with other orbitals, they give 
rise to a peak in the DOS. In contrast, 7r-bonding orbitals 
such as As p z and Pt d XZjVZ are located around the Fermi 
level with a wide band width. In particular, the van Hove 
singularity (vHS) coming from a saddle point near K just 
above Ep (inset of Fig. [5ja))exhibits two dimensional 
character since a 2D vHS gives rise to a singularity in the 
DOS while a 3D vHS gives only a discontinuity in slope. 
Raising Ep to vHS, which might be realized by electron 
doping via Sr layers, could potentially increase T c : As- 
suming rigid bands, we estimate that electron doping by 
roughly 25% will lift Ep to the vHS, enhancing the DOS 
at the Fermi level by 6%. Further, using a weak-coupling 
form for T c (T c = LUTne' 1 ^ ^), assuming that the 
pairing potential V does not change by electron doping, 
and assuming Debye temperature of Tp> = 229K^- we 
estimate that T c could be enhanced to as high as 3.8K 
by electron doping. 

Applying hydrostatic pressure could also raise the DOS 
and could, as often happens, enhance T c . Rh is a can- 
didate to supply a chemical pressure to SrPtAs since it 
has a similar electronegativity as Pt but a smaller ionic 
radius. However, the crystal structure is sensitive to the 
constituent atoms; for example, SrPtSb where As is re- 
placed by Sb has an AlBVtype structure, while YPtAs, 
where Sr is replaced by Y with more electrons, has a 
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hexagonal structure with four slightly puckered Pt As lay- 
ers in a unit cell£ 

Finally we consider the possibility of simple collinear 
magnetic solutions. The antiferromagnetic (AFM) phase, 
where the moments in a PtAs layer are aligned, by an- 
tiparallel to those in adjacent layers, is favored by 0.23 
meV (0.49 meV in GGA) per formula unit, than the non- 
magnetic phase; the ferromagnetic orientation converged 
to the non-magnetic solution. Magnetic moments are 
given in Table HU While the energy differences and cal- 
culated moments are too small to make definitive state- 
ments regarding magnetic phases in SrPtAs, the material 
appears to be near a magnetic instability. 



faces, i.e., a pocket and a cylinder, with 7r-character at 
the zone corner (around K-H). All the Fermi surfaces 
are hole-like which distinguishes this material from other 
pnictide superconductors. The transport properties are 
highly anisotropic between x, y— and z— directions. Rh 
is suggested for a positive pressure effect to increase T c . 
Furthermore, the van Hove singularity is shown in the 
DOS above Ep. Assuming rigidity of bands, we predict 
that Tc increases up 3.4 K with 25% doping, which may 
be achieved by chemical doping in place of the Sr atom. 
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TABLE II. Spin fi s and orbital \ii magnetic moments and 
total magnetic moment fit in unit of /is. 





Us 




Ht 


Pt 


0.014 


0.064 


0.078 


As 


0.006 


-0.027 


-0.021 
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